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We review the phenomenology of gravitino dark matter within supergravity 
framework. Gravitino can be dark matter if it is the lightest supcrsymmetric 
particle, which is stable if ij-parity is conserved. There are several distinct sce- 
narios depending on what the next to lightest supcrsymmetric particle (NLSP) 
is. We discuss the constraints and summarize the phenomenology of neutralino, 
stau, stop and sneutrino NLSPs. 
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1. Introduction 

The identity of dark matter has not yet been resolved at the time this 
article was written. One interesting possibility is to have a neutral lightest 
supersymmetric particle (LSP), which is stable via i?-parity conservation, 
as the dark matter. 1 Neutralino is a popular candidate for dark matter 
and sneutrino is another possibility a . Here, we focus our attention on yet 
another hypothesis within supergravity, i.e. gravitino as cold dark matter. 3 ' 4 
In supergravity models (i.e. models with gravity mediated supersym- 
metry breaking), 5 the gravitino mass is close to the other sparticle masses. 
However, it is not a priori whether the gravitino is lighter or heavier than the 
others. Note that this is different from gauge mediated models in which the 
gravitino mass can be naturally very light. 6 We assume supergravity models 
here, with supersymmetric masses of ~ 1 GeV — 1 TeV and the gravitino 
is the LSP. In this framework, the coupling between gravitino and mat- 



a Left-sneutrino is excluded by direct detection, 2 but right-sneutrino is still viable. 
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ter fields is very small, ~ 1/Mp\. Because of this, gravitino is practically 
undetectable (aside from its gravitational effect). Also, the next lightest 
supersymmetric particle (NLSP) could be long lived, with a lifetime of typ- 
ically 0(1 s) or longer. In this case, the NLSP decay affects the primordial 
light element abundances. 7 The phenomenology of this scenario depends 
largely on what the NLSP is. We will discuss below various possibilities, 
each with its own distinct phenomenology. 

2. Gravitino Dark Matter in Supergravity Models 

The biggest theoretical uncertainty in supersymmetric models arises from 
the fact that we do not know how supersymmetry, if it does exist, is broken 
in nature. Because of this, the values of the soft couplings are uncertain. We 
can take them as free parameters. However, because of the large number of 
parameters, we need to make some simplifying assumptions. Motivated by 
the Grand Unified Theory (GUT), the usual assumption is that parameters 
of the same type are unified at the GUT scale. Their values at weak scale 
are then derived by employing the renormalization group equation (RGE). 
The simplest model of this kind is the CMSSM (Constrained Minimal Su- 
persymmetric Standard Model) b , in which we have universal gaugino mass 
mi/2, universal sfermion mass too, and universal trilincar coupling Aq at 
the GUT scale. In addition, we have two parameters from the Higgs sector, 
i.e. the ratio of the two Higgs vevs tan/3 = (Hi)/(H 2 ) 7 and the sign of /i 
where /x is the Higgs mixing parameter in the superpotcntial. Note however 
that in GUT theories, e.g. SU(5) or 5*0(10), the Higgs fields are contained 
in different multiplets as compared to the matter multiplets. This moti- 
vates a generalization of the CMSSM, in which the Higgs soft masses 7711,2 
are not necessarily equal to to at the GUT scale. 8 Furthermore, we can 
trade toi j2 with /x and the CP-odd Higgs mass 771,4 as our free parameters 
through the electroweak symmetry breaking condition. The resulting model 
is called Non-Universal Higgs Masses (NUHM) model. 9 Thus the NUHM 
parameters are 777,1/2, ™o, Ag, tan/3, \x and rriA- 

In the usual scenario with neutralino LSP, we implicitly assume that 
gravitino is sufficiently heavy such that it decouples from the low energy 
theory. In the gravitino dark matter (GDM) scenario, on contrary, we as- 
sume that the gravitino mass = 777,3/2 is sufficiently small such that the 
gravitino is the LSP. For our purposes, we can take 777,3/2 as another free 
parameter. Within CMSSM, with gravitino LSP, there are three possible 



b Also known as mSUGRA, depending on your preference. 
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NLSP, i.e. neutralino, stau and stop particles. 3 ' 10 For NUHM, in addition, 
we can have selectron or sneutrino as the NLSP. 11 Of course for a more 
general model of MSSM we can have more possibilities. 

3. Phenomenological Constraints 

3.1. Dark matter relic density constraint 

Being a very weakly interacting particle, gravitino decoupled very quickly 
from the thermal plasma in the early universe. This leads to a concern that 
the gravitino could be over-abundance. However, inflation can solve this 
problem. 12 In inflationary models, the early gravitino density together with 
other densities are diluted by the inflation. Gravitino is then reproduced by 
reheating after the inflation, 13 although with a smaller yield that can still 
satisfy the relic density constraint c . 

Gravitino relic density consists of two parts, the thermal relic Q~ which 
is produced by reheating, and the non-thermal relic f2 ~ T coming from the 
decay of the NLSP. 

n d h 2 = n T d h 2 + n^ T h 2 (i) 

The thermal relic is related to the reheating temperature T R through the 
following relation 15 




where m g is the gluino mass. We can see that for — 100 GeV and 
m„ = 1 TcV, to get n T ~h 2 < 0.1 we need T R < 10 10 GeV. The value of 
T R depends on the inflation model. For our purpose, we take T R as a free 
parameter. For one to one decays of NLSP to gravitino, which is generally 
the case, the gravitino non-thermal relic can be written as: 

nfh 2 = -^jw^ ( 3 ) 

where ^nlsp^ 2 is the NLSP density before the decay. Due to the long 
lifetime, the NLSP density is frozen out long before its decay, and this can 
be calculated by the usual method of solving the Boltzmann equation in the 
expanding universe. Note that even if the NLSP density is larger than the 
WMAP value, we might still satisfy the dark matter relic density constraint 



c This can still impose a strong constraint on the inflationary theories. 14 However, this 
topic is beyond the scope of this article. 
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because of the rescaling by the mass ratio mg/mNLSP- The NLSP density 
can also be written in term of the yield, Inlsp = ^nlsp/ s 7 where n is the 
number density and s is the entropy density. This is related to ^nlsp^ 2 by 

1nlspM NLS p = ft NLSP /i 2 x (3.65 x 1(T 9 GeV) (4) 

The total relic density of the gravitino must not exceed the upper limit 
of the dark matter relic density range as suggested by WMAP: 16 

Cf P h 2 ~ 0.113 ± 0.004 (5) 

Taking 2cr, this means that Q^h 2 < 0.121. The percentage of the thermal 
versus non-thermal relic density depends on the strength of the NLSP in- 
teractions which determine il^LSph 2 - If we can measure these interactions 
at colliders we can deduce the reheating temperature. 17 

3.2. The BBN constraints 

Big bang nucleosynthesis (BBN) is often cited as the greatest success of 
the big bang theory. By using simple assumptions that the early universe 
is in thermal equilibrium and expanding one can calculate the primordial 
light clement abundances (using standard nuclear cross sections) and gets 
results which agree very well with the observations. If there is a metastable 
particle that decays during or after the BBN era, the light element abun- 
dances can be altered by the participation of the energetic decay products 
in the nucleosynthesis processes. Thus, BBN provides a stringent constraint 
for gravitino dark matter. On the other hand, the prediction of the stan- 
dard BBN (sBBN) is not perfectly in agreement with the observational 
data. There seems to be discrepancy between the observed lithium abun- 
dances and the predicted values as shown in Table 1. This is known as 



Table 1. Comparison of lithium abundances from standard model 
prediction and observations. 





7 Li/H 


6 Li/ 7 Li 


Observation 


(1 - 2) x 10 -10 


~ 0.01 - 0.15 


Standard BBN (with CMB) 


~4x 10" 10 


< io- 4 



the lithium problem. The discrepancy on 6 Li is particularly difficult to be 
solved within the standard theory. There is no known astrophysical process 
that can produce large amount of 6 Li. Moreover 6 Li is fragile. Therefore 
we should expect less rather than more 6 Li compared to the prediction. 
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The lithium problem could be an indication of a new physics beyond the 
standard model. There are two proposed solution to this problem through 
a hypothesized metastable particle. The first one is through catalytic ef- 
fect. 18 ' 19 The process d+ 4 He — ► 6 Li +7 is suppressed by parity. If there is 
a massive negatively charged particle X~ that is bound to 4 He by Coulomb 
interaction it can absorb the emitted photon, hence the process is no longer 
parity suppressed. Simultaneously, the X~ particle is freed from the bound 
by the energy released and can subsequently be attached to another 4 He, 
thus acting as a catalysis for 6 Li production. This catalytic effect can also 
effect other light element abundances such as beryllium. 20 

Another proposed solution to the lithium problem is through hadronic 
decay of a metastable particle. 21 The decay produces energetic n, p and 
also T, 3 He (through spallation of 4 He), which then interact with the am- 
bient nuclei, e.g. n + p — >D, T + 3 He -^ 6 Li (producing more 6 Li), and 
7 Be(n,p) 7 Li(p, 4 He) 4 He (reducing 7 Li). Note that deuterium is also pro- 
duced, hence put some constraints on this scenario. 

3.3. Astrophysical constraints 

If the NLSP decays at a time later than the BBN, the photons produced by 
the decay might not be fully thcrmalised by the time of recombination. This 
can cause distortion on the cosmic microwave background radiation black 
body spectrum. 22 The size of the distortion depends on the amount of the 
energy injected into photons. This is represented by a chemical potential \x 
for photon. CMB spectrum measurement by the FIRAS instrument onboard 
of COBE satellite sets an upper limit on fi: 23 

\n\ < 9 x 1CT 5 (6) 

This limit is only important for lifetime > 10 6 s since photons produced 
earlier should have enough time to thermalize before recombination. 

Gravitinos from the NLSP decay at a late time have larger velocities 
compared to the primordial gravitino. This leads to a longer free-streaming 
length, smoothing out small scale density perturbation. If the dark matter 
relic density is dominated by non-thermal relic, the structure formation 
would be affected. This scenario is proposed as a solution to the small scale 
problem. 24 

3.4. Collider constraints 

At colliders, heavier supersymmetric particles can still be produced pro- 
vided that there is enough energy in the collisions. These sparticle would 
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quickly decay, cascading down to the NLSP. Due to its long lifetime the 
NLSP itself would escape from the detector before eventually decays to 
gravitino, hence appear as a stable particle with respect to the detectors. 
There have been searches for stable massive particles (SMP) at colliders. 25 

A particularly interesting signal would be produced if the NLSP is 
electromagnetically charged. In this case the NLSP should traverse the 
calorimeter and subsequently be detected by the muon detector. The first 
obvious step of the data analysis is of course to discover this charged NLSP. 
The CDF collaboration, based on 1.0 fb _1 of data at y/s = 1.96 TeV, sets 
a lower bound on (meta)stable stop particle at 249 GeV; 26 while the DO 
collaboration, using 1.1 fb _1 of data, sets upper limits for stable stau pair 
production cross section from 0.31 pb to 0.04 pb for stau masses between 
60 GeV and 300 GeV 27 and lower mass limits of 206 GeV and 171 GeV for 
pair produced stable charged gauginos and higgsinos respectively. 

However, not all possible NLSP are charged. If neutralino or sneutrino 
is the NLSP, they would not be detected. Only neutralino with a very short 
lifetime (< few ns) can be detected through its decay product. The CDF sets 
a lower limit on the neutralino mass at 101 GeV for lifetime 5 ns. 28 Similar 
to the familiar case with neutralino LSP, there are also various signatures 
from the cascade decays. The same methods of analysis can be applied 
to the case with stable neutral NLSP. For long-lived neutralino NLSP the 
signatures would be indistinguishable from that of neutralino LSP scenario. 
For sneutrino NLSP, however, the signatures would in general be different. 29 

4. Phenomenology of GDM with Various NLSP 

In this section we look at each scenario of neutralino, stau, stop and sneu- 
trino NLSP. We do not include chargino NLSP 30 here. 

4.1. Neutralino NLSP 

For neutralino mass of 1 TeV and gravitino mass of 1 GeV the neutralino 
lifetime is about 0(1) s. The lifetime is longer for a smaller mass gap. 
Thus the neutralino in this scenario would escape the collider detectors 
and trigger large missing energy signatures. Assuming that all primordial 
neutralinos has eventually decayed to gravitino, only gravitino is floating 
around today. Therefore WIMP direct detection experiments would not 
see any signal. There would be no indirect astrophysical signal from dark 
matter annihilation in the halo either. In this case it would be difficult to 
proof the identity of the dark matter, i.e. whether it is gravitino, axino or 
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whether something else. We will also need to check whether i?-parity is 
really conserved. 

The neutralino NLSP in the CMSSM is much constrained by the BBN, 
especially when the neutralino lifetime is > 10 4 s where the electromagnetic 
shower effect on the light element abundances becomes important. The 
main reason is because neutralino has relatively large freeze-out density for 
most of the parameter space. 

4.2. Stau NLSP 

If produced, a stau NLSP would be seen as a massive stable charged particle 
at colliders. It would leave a clean track in the inner detector and then reach 
the muon detector, hence it would look like a slow/heavy muon. Because 
of its electromagnetic charge, the stau can be slowed down by making it 
go through a bulky medium. Thus it can be trapped and stored until it 
decays. 31 ' 32 In this way one can hope to measure its lifetime. 

Within the CMSSM, stau NLSP has the largest allowed region of pa- 
rameter space, hence thought as the natural candidate. Stau NLSP would 
yield catalytic effect on BBN. This has attract much attention and many 
papers are devoted in the study of this topic, in particular regarding the 
lithium problem solution. 33 

4.3. Stop NLSP 

A long lived stop would hadronizc once it is produced. By taking anal- 
ogy with heavy quark hadrons one can deduce the lightest hadron states 
and their lifetimes. The light stop sbaryons are A~ = t\ud (which is the 
lightest), = ti(uu,ud,dd) (which decays through strong interac- 

tion), and = tis(u,d) (which decays semileptonically with lifetime 

t < 10~ 2 s). The light stop mesinos arc T° = i\u (which is the light- 
est), T+ = t\d (with lifetime r ~ 1.2 s), and T s = t\s (with lifetime 
t~2x 10 -6 s). The antistop would hadronize into the corresponding an- 
tisbaryons and antimesinos. In the early universe, being the lighter one, 
the neutral T° is more abundance than the charged A~. This reduces the 
catalytic effect on BBN. Moreover, due to its strong interaction, the freeze 
out density of stop is generally small. Therefore this scenario can generally 
satisfy the BBN constraint. 10 

On the other hand, it was shown 34 that stop NLSP scenario is fit to solve 
the lithium problem through hadronic decay. Note that further annihilation 
of stop occurs after the hadronization, with annihilation rate L ann = (av)rif 
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where a ~ -Rhad an( ^ v — V^/ m t- The final stop abundance before its 
decay can be written as 

^ = 0.87xlO-Gev(A)- 2 (_^_)- 1/2 

x ( r ^D Y 3/2 ( m y /2 a) 

Vl50McVy Vl0 2 GcV/ V ; 

It was found that, with f a = 0.1, the lithium problem solution prefers 
m t - = 400 - 600 GeV and m 3/2 = 2-10 GeV. 

4.4. Sneutrino NLSP 

Since sneutrino does not interact strongly nor electromagnctically the effect 
of sneutrino decays on BBN can be guessed to be small, but nonzero. The 
BBN effect comes through energy transfer (clastic/inelastic) from energetic 
neutrino (produced by the decay v — > G + v) to the background particles, 
and from 3/4-body decay modes 35 [y — » G + v + (7, Z), v — > G + i + W) 
and 4-body {v -> G + v + / + /,/> -» G + £ + / + /')• Although the 3 and 
4 body decay branching ratios are small, they can still be important since 
they produced particles that can directly involve in the nucleosynthesis 
processes. For Mp <~ O(100 GeV), the BBN constraint can be satisfied if 
the sneutrino density before the decay is 

YpMs < 0(lO~ n ) GeV for B h = 10~ 3 (8) 

YvMp < O(10- s ) GeV for B h = 10" 6 (9) 

where Bh is the hadronic branching ratio. The sneutrino NLSP and grav- 
itino LSP scenario was explored within NUHM models, and it was found 
that there are large regions of parameter space still allowed. 11 

At colliders, similar to the neutralino (N)LSP case, sneutrino NLSP 
would yield a missing energy signature. We can study this scenario through 
cascade decays of heavier supersymmctric particles. In general, the signa- 
tures are different from those in neutralino case. Signatures that are thought 
to be the best for neutralino LSP might not be suitable for this sneutrino 
NLSP case. A preliminary study of collider phenomenology with sneutrino 
NLSP has been done in Ref. 29. However, a more detail study might reveal 
more information. 

5. Concluding Remarks 

Gravitino is a feasible and interesting candidate for dark matter. There are 
many possible phenomenology with gravitino dark matter depending on the 
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choice for the NLSP, which we have summarized in this article. Future and 
upcoming collider experiments, such as the LHC, might be able to unveil 
some hints on the identity of dark matter. Progresses in direct and indirect 
detection experiments are also looked promising. The next few years would 
be an interesting time to find out more about dark matter, and whether 
gravitino can still stand up as a candidate for dark matter. 
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